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Preface 

Even  a  casual  reading  of  newspaper  or  magazine  articles 
suggests  that  ours  is  a  world  which  is  vitally  concerned  with 
both  the  production  and  conservation  of  energy.   People  every- 
where are  searching  for  new,  cheaper,  and  environmentally 
safe  energy  resources,  while  others  are  trying  to  find  more 
efficient  ways  of  using  the  energy  we  already  have. 

This  work  is  consistent  with  those  efforts  in  that  it 
seeks  to  further  develop  a  technology  which,  properly  applied, 
could  contribute  to  the  solution  of  both  aspects  of  the  energy 
problem.   To  be  specific,  we  are  here  concerned  with  developing 
a  procedure  for  the  numerical  design  of  transonic   turbine 
and  compressor  blade  profiles   in  two  dimensions.    We  say 
that  a  flow  is  transonic  if  it  attains  speeds  both  above  and 
below  that  of  sound.   Such  a  procedure  might  be  used  to  design 
more  efficient  axial  flow  compressors  for  use  in  the  produc- 
tion of  enriched  uranium  at  gaseous  diffusion  plants,  as  well 
as  in  the  construction  of  lighter,  more  efficient  airplane 
engines  for  better  fuel  consumption.   Improving  the  design  of 
high  speed  turbomachinery  could  be  of  practical  benefit  to  a 
great  many  people. 

The  techniques  used  here  are  a  direct  result  of  work 
on  the  supercritical  wing.  First  designed  by  Whitcomb  and 
his  associates  at  the  Langley  Research  Center  in  Virginia, 
such  airfoils  have  undergone  wind-tunnel  and  flight  testing. 
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the  results  showing  that  it  is  indeed  possible  to  achieve 
virtually  shock-free  flow  without  boundary  layer  separation. 
As  a  consequence  of  these  experimental  successes,  attempts 
at  providing  a  purely  theoretical  basis  for  the  design  of 
supercritical  wings  have  become  more  frequent  and  increas- 
ingly successful.       In  fact,  wing  sections  designed  by 
the  method  used  here  have  been  tested  also  [17],  and  have 
shown  the  theory  to  be  quite  functional . 
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Introduction 

1.    The  Physical  Problem 

In  this  section,  we   use  the  gas  turbine,  with  particular 
attention  to  its  use  in  jet  propulsion,  to  provide  a  physical 
context  for  the  mathematical  problem  under  study. 

In  general,  the  gas  turbine  exists  in  a  variety  of  forms 
and  has  many  applications,  though  its  operating  principle  is 
quite  simple.   In  its  basic  form,  the  compressor  section  raises 
the  pressure  and  temperature  of  the  gas,  the  combustion  section 
adds  fuel  and  oxidizes  it,  thus  raising  the  temperature  sharply, 
and  the  turbine  section  fans  it  out,  lowering  the  pressure 
considerably  and  the  temperature  slightly.   The  useful  work 
of  the  system  is,  of  course,  delivered  to  the  turbine  wheel. 
Note  that  because  of  the  work  done  in  compressing  the  gas,  the 
compressor  is  usually  driven  by  the  turbine. 

When  used  as  a  power  plant,  the  turbine  wheel  delivers 
only  a  part  of  its  energy  to  the  compressor  (when  there  i£  a 
compressor),  the  remainder  being  used  elsewhere.   In  a  jet 
propulsion  system,  however,  all  of  the  turbine  wheel  energy 
goes  to  the  compressor,  the  momentum  of  the  craft  being  a 
result  of  Newton's  third  law  and  the  rearward  discharge  of 
gases.   A  much  simplified  diagram  of  such  a  machine  appears 
in  Figure  1.   Observe  that  though  there  are  several  types  of 
compressors,  we  will  at  all  times  be   restricted  to  axial 
flow  compressors,  so  nam.ed  because  the  gas  flows  along  the 
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axis  (see  Figure  2) . 

If  we  now  look  at  a  blade  row  of  such  a  structure, 
or  a  turbine  wheel   (see  Figure  3) ,   we  see  that, 
though     quite  different,  each  consists  of  a  pattern  of 
blades  mounted  on  a  right  circular  cylinder.   Further,  we 
ignore  the  radial  components  of  the  gas  velocity,  i.e.  we 
in  essence  assume  that  the  flow  streamlines  lie  on  the 
cylinder's  surface,  as  in  Figure  4.   Finally,  we  develop 
this  cylindrical  surface  into  a  plane,  obtaining  the  two 
dimensional  cascade  in  Figure  5. 

These  are  the  structures  we  want  to  design.   More 
precisely,  we  wish  to  find  smooth,  transonic,  steady-state 
solutions  to  the  equations  for  the  irrotational  flow  of  an 
inviscid  fluid  through  a  two  dimensional  cascade. 

It  should  be  pointed  out  that  industry  has  found 
cascades  to  be  quite  useful  in  blade  design  work,  even  though 
they  are  only  two  dimensional  idealizations  of  the   physical 
fact. 
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2.    The  Mathematical  Problem 

In  mathematical  terms ,  our  problem  reduces  to  finding 
analytic  solutions  to  a  system  of  partial  differential  equa- 
tions for  flow  about  a  body.   Although  we  will  often  be  concerned 
with  a  single  blade,  the  periodicity  of  the  solution  will 
result  in  a  cascade.   There  are  several  difficulties  in  such 
a  problem: 

(1)  The  equations  are  nonlinear, 

(2)  Analytic  solutions  may  not  always  exist,  and 

(3)  The  Equations  are  of  mixed  type. 

We  overcome  the  first  difficulty  by  use  of  the  hodograph 
transformation.   First  used  by  Chaplygin  in  his  work  on  gas 
jets  [  6  ] ,  this  well-known  procedure  interchanges  the  roles 
played  by  the  independent  variables  (the  planar  coordinates 
X  and  y)  and  the  dependent  variables  (the  velocity  components 
u  and  v) .   The  resulting  equations  are  linear,  though  the 
boundary  data  becomes  somewhat  complicated. 

The  second  difficulty,  as  well  as  the  third,  are  results 
of  the  transonic  nature  of  the  flow  and  can  be  presented  in  the 
following  way.   Suppose  we  wish  to  design  a  turbine  blade. 
One  procedure  is  to  input  a  particular  blade  design  with  the 
relevant  informatior\  e.g.   the  geonetry  of  the  blade  surface, 
and  compute  the  flow.   This  is  called  the  direct  problem, 
since  we  specify  the  blade  and  compute  the  flow.   It  has  been 
shown  by  Morawetz  [25 1,  however,  that  the  problem  of  finding 
a  shock-free  transonic  flow  about  an  arbitrary  body  is  not 
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well-posed,  i.e.  a  smooth  solution  may  not  exist.   We  can  get 
around  the  issue  by  instead  solving  the  inverse  problem,  i.e. 
we  choose  data  which  specify  properties  of  the  flow,  and  then 
compute  the  blade  profile.   Thus  the  solution  will  be   analytic, 
and  we  can  then  manipulate  the  data  until  we  get  a  good  blade. 

The  remaining  difficulty  is,  essentially,  that  of 
finding  an  applicable  method  for  solving  the  equations.   Like 
Tricomi ' s  equation,  the  equations  of  transonic  flow  are  of 
mixed  type  (hyperbolic  and  elliptic)  and  thus  require  special 
algorithms  for  their  solution.   The  procedure  chosen  here  is 
the  method  of  characteristics.   This  technique  is  usually 
restricted  to  hyperbolic  systems,  but  we  have  used  it  every- 
where, after  extending  all  variables  into  the  complex  domain 
by  analytic  continuation.   The  resulting  equations  are  then 
cast  in  the  form  of  a  complex  characteristic  initial  value 
problem  and  solved  in  complex  four-space.  Since  we  are  using 
the  hodograph  equations,  which  are  linear,  we  can,  by  the 
principle  of  superposition,  break  up  our  expected  solution 
into  two  parts,  a  regular  part  and  a  singular  part,  the  form 
of  the  latter  being  dictated  by  the  cascade  geometry. 

The  study  of  transonic  problems  received  a  great  impetus 
as  a  result  of  experimental  successes  in  achieving  virtually 
shock-free  transonic  flow  in  the  1960 's  [19],  [27],  [29]- 
The  particular  technique  of  treating  these  problems  as  complex 
characteristic  initial  value  problems,  first  used  by  Garabedian 
and  Lieberstein  [11 ] ,  has  been  carried  forward  to  the  design 
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of  airfoils  [  2] ,  [  1] .   The  initial   application 
to   compressors   is   documented   by   Korn   in   [21] . 
It  should  be  mentioned  that  weak  shocks  appear  under  off- 
design  conditions,  though  the  losses  remain  low. 

Nevertheless,  certain  difficulties  have  arisen.  It 
has  been  found  difficult  to  generate  blades  of  large  thick- 
ness and  high  camber.  The  contribution  of  this  work  has  been 
to  design  such  blades.  This  was  done  by  basing  the  form  of 
the  characteristic  mapping  on  the  complex  velocity  of  a  flow 
about  a  Joukowski  blade  mapped  onto  the  interior  of  a  circle, 
instead  of  the  complex  velocity  of  a  circular  cylinder  as  had 
been  done  before. 

The  plan   of  this  paper  is  as  follows.   In  Chapter  1 
we  use  basic  physical  assumptions  to  derive  the  characteristic 
equations  and  extend  them  into  the  complex  domain,  the  result 
being  a  system  of  equations  which  apply  to  any  transonic  flow. 
In  Chapter  2,  we  present  a  singular  solution  which  is  based 
on  the  geometry  of  a  cascade  and  describe  the  selection  of 
appropriate  initial  data,  the  result  being  a  threefold  complex 
initial  value  problem  which  applies  to  a  transonic  cascade. 
Chapter  3  contains  a  descriptive  account  of  the  numerical 
method  of  solution  (integration  paths,  finite  difference 
scheires,  graphic     output,  etc.),  emphasizing  the  basic 
strategy  of  the  computer  program,     as   well   as 
the  results  of  the  original  program,  viz.,  that  we  were  able 
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to  generate  a  compressor  blade  but  not  a  turbine  blade.   Last, 
in  Chapter  4  we  introduce  the  Joukowski  transformation  and 
present  the  resulting  turbine  blade. 

There  are,  of  course,  many  people  who,  at  one  point  or 
another,  helped  me  in  this  work.   In   particular,  I  would 
like  to  acknowledge  my  tremendous  debt  to  Paul  R.  Garabedian, 
for  his  many  suggestions,   tireless  support,  and  enduring 
patience.   He  suggested  this  problem,  and  it  goes  without 
saying  that  the  work  could  not  have  been  done  without  him. 
I  am  also  grateful  to  David  Korn  for  his  aid  and  encouragement, 
and  for  allowing  me  to  use  his  programs.   Thanks  are  also 
due    Connie  Engle   for   typing  the  manuscript. 
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1.  TKE  PARTIAT.  DIFFERENTIAL  EQUATIONS  OF  TRANSONIC  FLOW 

In  this  chapter  we  state  the  basic  equations  of  the 
flow  and  introduce  the  hodograph  transformation.   The 
resulting  equations,  which  are  of  mixed  type,  are  then  cast 
in  the  form  of  a  characteristic  initial  value  problem  in 
the  complex  domain. 

1.     The  Flow  Equations 

The  mathematical  model  for  the  inviscid  flow  of  a  com- 
pressible fluid  is  based  on  the  following  assumptions  ( [  7  ] , 
p.  3-25],  as  expressed  in  the  accompanying  equations. 
Note  that  we  will  use  Euler ' s  representation   at  all  times. 

(1)    The  conservation  of  mass: 


p^  +  V- (pq)  =  0  , 


where  p,  the  density,  and  q,  the  fluid  velocity,  are  functions 
of  space  and  time,  p   is  the  partial  derivative  of  p  with 
respect  to  t,  and 

V  =  i  I-  +  j  I-  +k  I-  . 
9x    -^  3y     3z 

(2)    The  conservation  of  linear  momentum: 


P  dt 


g  +  VP  =  0  , 
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where  P  represents  the  pressure,  and 


dt    3t      3x      3y      3z 


is  the  material  derivative.   Here  and  in  the  following  we 
will  let  u,  V  and  w  respectively  denote  the  x,  y  and  z 
components  of  the  velocity  vector  q. 

(3)    The  equation  of  state  of  a  polytropic  gas: 


P  =  A(S)p^  , 


where  S  is  the  specific  entropy,  A  is  a  known  function  of  S, 
and  Y  is  the  adiabatic  exponent. 

(4)    The  conservation  of  energy: 

^=  0 
dt     • 

This  last  equation  is  somewhat  unnecessary  for  our  particular 
purposes  since  we  assume  that  the  flow  is  isentropic,  i.e.  S 
is  constant  throughout  the  medium. 

If,  in  addition,  we  restrict  ourselves  to  the  steady 
flow  of  a  polytropic  gas  in  two  dimensions,  the  remaining 
equations  reduce  to 

(1.1)  (pu)^  +  (pv)   =  0 

(1.2a)  p(uu   +  vu  )  +  P   =0 

X     y     X 

(1.2b)  p  (uv   +  w  )  +  P   =0 

X     y     y 
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and 


(1.3)  P  =  Ap^, 


where  A  is  now  a  known  constant.   Observe  that  we  can  differ- 
entiate (1.3)  to  get 


(1.4)  c'  =  J-  =  YAp 


2  _  dP y^^y-i 

dp 


which  defines  the  local  speed  of  sound  c.   Consequently  we 

can  use  both  (1.3)  and  (1.4)  to  eliminate  p  and  P  from  (1.1)  and 

(1.2)^  and  then  manipulate  the  resulting  equations  to  get 

(1.5)  (c^-u^)u^  -  uv(u  +v  )  +  (c^-v^)v   =  0  . 

X       y   X  y 

Finally,  we  assume  that  the  flow  is  irrotational,  i.e. 

(1.6)  "y  "  ^x  "  °  ' 

which,   together  with  (1.2),  implies  that 

2    2c2 

(1.7)  q   +  TTi  =  TTI  ' 

where  c-  is  the  sound  speed  at  q  =  0.   Equation  (1.7)  is, 
of  course,  a  form  of  Bernoulli's  law. 

In  summary,  we  remark  that  the  system  of  equations  (1.5) , 
(1.6)   corresponds  to  the  steady,  isentropic,  irrotational 
flow  of  a  polytropic  gas  in  two  dimensions.   We  wish  to  solve 
this   system  for  u  and  v,  where  we  consider  c  to  be  a  known 
function  of  q  given  by  (1.7). 
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2.    The  Hodograph  Equations 

Clearly,  equation  (1.5)  is  nonlinear,  since  the  coeffi- 
cients multiplying  the  partials  of  u  and  v  are  themselves 
functions  of  u  and  v.   One  may,  however,  by  use  of  implicit 
function  theorems  or,  equivalently ,  a  Legendre  transformation 
{[7],  pp.  248-252),  rewrite  (1.5)  and  (1.6)  as 

(1.8)  (c^-u^)y^  +  uv(x^+y^)  +  (c^-v^)x^  =  0 
and 

(1.9)  X   -  y  =  0  . 

V   -'u 

Of  course,  we  have  assum.ed  that  the  Jacobian 

3(x,y)  ^ 

in  the  region  of  interest. 

The  important  point  here  is  that  by  considering  x  and  y  to 
^^  functions  of  the  now  independent  variables  u  and  v  we 
obtain  a  linear  system  of  first  order  partial  differential 
equations.   Equations  (1.8)  and  (1.9)  are  called  the 
hodograph  equations . 
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3 .     The  Characteristic  Equations 

We  first  rewrite  the  systen  (1.8),  (1.9)  as 


(1.10) 


SZ   +  TZ   =  0  , 
V      u 


where  we  have  set 


S  = 


T  = 


2   2 

c  -u 


uv 

-1 


and 


Z  = 


y 

X 


uv 


2   2 

C  -V 


We  wish  to  solve  this  system  by  the  method  of  characteristics. 
To  do  that,  we  first  find  equations  in  the  uv-plane  for  the 
two-pareimeter  family  of  loci  called  characteristics,  or 
characteristic  curves,  the  two  parameters  being  called 
characteristic  coordinates.   We  then  construct  equations 
for  the  rem.aining  variables,  y  and  x,  in  terms  of  the  charac- 
teristic coordinates,  and  solve  the  entire  system  numerically 
([ 7  ] ,  pp.  40-45) . 

Toward  that  end,  we  consider  the  smooth  curve  C  contained 
in  the  uv-plane  and  given  by  (u(s),v(s))  where  s  is  a  real 
parameter.   Then,  using  the  relations 


dZ  =  y  ^  +  y  ^  , 
ds    v  ds   ■'u  ds 


and 


-11- 


dx  _    dv      du 
ds    V  ds     u  ds  ' 


we  may  rewrite  (1.10)  in  the  form 


(1.11) 
where 


AX  =  B  , 


A  = 


2      2 
c   -u 

uv 

uv 

2      2 

C     -V 

0 

1 

-1 

0 

dv 
ds 

0 

du 
ds 

0 

0 

dv 
ds 

0 

du 
ds 

X  = 


■   ^v    ' 

^v 

^u 

X 

u 

and 


0 

d^ 
ds 

dx 
ds 

We  recall  that  in  two  dimensions,  however,  a  characteristic 

is  a  curve  along  which  one  cannot  solve  the  given  equations 

for  the  normal  derivatives  of  the  dependent  variables  (in 

this  case  y  and  x) ,  even  though  the  values  of  these  functions 
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are  known  on  that  curve.   Nevertheless,  if  we  know  y  and  x 
on  C,  we  can  differentiate  them  to  get  dy/ds  and  dx/ds  there, 
which  means  that  we  can  determine  the  matrices  A  and  B. 
Further,  if  A  is  nonsingular,  we  can  invert  it  and  get  a 
unique  X,  which  means  that  we  will  have  determined  the 
derivatives  of  y  and  x  in  all  directions,  including  the 
normal.   Hence,  if  C  is  a  characteristic,  A  is  singular,  i.e. 

(1.12)  det  A  =  0  . 

Introducing  X  =  dv/du  ,  we  assume  du/ds  /  0  and  set 


<l-l^'  37=^  if 


Then,  carrying  out  the  computation  of  (1.12)  we  get 

2 

det  A  =  (J^)  [X^(c^-v^)  -  2uvX  +  c^-u^]  =  0  , 
as 


which  implies  that 


(1.14)         X^(c^-v^)  -  2uvX  +  c^-u^  =  0 


It  is  this  equation,  in  practice,  which  is  used  to  determine 
the  characteristic  curves,  since  we  do  not  know  their  equa- 
tions beforehand. 

Solving  (1.14)  for  X,  we  get  two  solutions. 


4.   /  2     2 
uv  ±  c/a   -  c 
(1.15)  X+  =  — ^ f—   , 

c   -  V 
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corresponding  to  two  families  of  characteristic  curves  which 
we  shall  hereafter  refer  to  as  C   and  C_.   We  call  X^(u,v) 
the  characteristic  directions,  since  they  are  the  slopes  of 

the  characteristics  through  each  point  of  the  uv-plane. 

2     2 
Clearly  they  are  real  and  distinct  whenever  q   >  c  ,  i.e. 

whenever  the  flow  is  supersonic.   Thus  we  have  determined  a 

two  parameter  family  of  curves  in  the  uv-plane  along  which 

we  wish  to  solve  for  y  and  x. 

Suppose  now,  that  we  can  find  two  parameters,  a  and  B, 

such  that 

B(u,v)  =  constant 
on  each  curve  in  C   ,  and 

a(u,v)  =  constant 

on  each  curve  in  C_  ,   where 

3(a,e) 


3(u,v) 


i^      0  . 


We  then  call  a  and  B  characteristic  coordinates,  since  their 
values  at  a  point  in  the  uv-plane  indicate  which  character- 
istics intersect  there.   Thus  we  may  rewrite  (1.13)  as 

(1.16a)  v^  =  X^u^ 

on  C_|^  ,  and 
(1.16b)  Vg  =  X_Ug 

on  C   ,  where  we  have  used  partial  derivatives  since  u  and  v 

(as  well  as  y  and  x)  are  assumed  to  be  functions  of  both  a  and  B. 
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Observe  then,  that  in  the  case  of  C   ,  (1.11)  leads  to 
the  augmented  matrix 


[A,B]  = 


2      2 
c    -u 

uv 

uv 

C^-v 

2 

0 

0 

1 

-1 

0 

0 

+   a 

0 

u 
a 

0 

^a 

0 

X    u 

+   a 

0 

"a 

"a 

which,  by  use  of  elementary  row  operations,  can  be  reduced  to 


0 

0 

0 

0 

^a+^-^a    ' 

0 

1 

-1 

0 

0 

+ 

0 

1 

0 

^a/^a 

0 

0 

K 

1 

'^a/^a 

The  last  three  rows  yield  no  new  information,  but  from  the 
first,  we  see  that  we  must  have 


(1.17a) 


X   +  A  y   =0 
a    -■'a 


on  every  curve  of  C   ,  or  our  system  will  be  inconsistent. 
In  a  similar  way,  we  arrive  at 


(1.17b) 


Xg  +  X^y  =  0 


for  every  curve  in  C_.   Thus,  equations  (1.16),  (1.17)  form 
a  system  of  four  characteristic  equations  in  four  unknowns 
(u,v,y  and  x) ,  in  terms  of  the  as  yet  undetermined  character- 
istic coordinates  a  and  B. 
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Explicit  expressions  for  a  and  g  can  be  found,  however; 
in  fact.  Sears  has  shown  [28 ]  that 

(1.18a)  a  =  f (q)  -  6 

{1.18b)  e  =  f(q)  +  e  , 

is  such  a  pair  of  parameters ,  where  of  course 

/T~   2 
q  =  /u   +  V 


-Irvi 


e  =  tan'-'H  , 


and  where 


(1.19)   f (q)  =  ^  sin"^ 


2 

(Y-1)  ^  -  Y 
c 


^  1   .  -1 
+  2   sin 


c* 


(Y+1)  -^  -  y 


* 


Here,  c^  is  the  critical  speed,  and 


^'/U  ■ 


Note  that  characteristic  coordiates  are  certainly  not  unique, 
since  any  continuously  dif f erentiable  monotone  function  of 
a  characteristic  coordinate  is  itself  a  characteristic  coordi- 
nate.  Our  particular  choice  of  characteristic  coordinates  (see 
Section  2.2)  will  be  based  on  Sears',  though  they  will  be,  in 
general,  complex,  and  have  additional  properties  as  consequences 
of  that  fact. 

In  passing,  we  remark  that  we  would  have  obtained  equa- 
tions (1.16),  (1.17)  even  without  the  hodograph  transformation. 
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with  -1/X+  as  the  characteristic  directions  in  the  xy-plane. 
Our  approach  differs  in  that  we  will  solve  (1.16)  for  u  and 
V  in  terms  of  a  and  3,  and  then  prescribe  initial  data  for 
y  and  x  in  the  uv-plane.   The  advantage  gained  is  that  equa- 
tions (1.17)  are  linear,  since  A^  are  functions  of  u  and  v 
alone,  and  we  can  use  the  principle  of  superposition. 
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4 .     Complex  Extension 

The  distinguishing  characteristic  of  a  transonic  flow  is 
that  it  contains  regions  in  which  the  flow  speed  is  above  the 
local  speed  of  sound  as  well  as  regions  in  which  it  is  below. 
In  the  former  case,  i.e.  when  q  >  c,  the  system  is  hyperbolic 
and,  since  the  characteristic  directions  X^  are  real  and 
distinct,  we  can  use  the  method  of  characteristics.       When 
the  flow  is  subsonic,  however,  the  characteristic  directions 
take  on  complex  values  and  the  method  breaks  down.   The  system 
is  elliptic  in  such  a  region,  and  the  characteristic  curves 
are  not  defined. 

We  resolve  this  issue  by  extending  all  variables  into 
the  complex  domain  and  solving  the  original  two  dimensional 
problem  in  complex  four-space.   That  this  can  be  done  is  a 
consequence  of  the  fact  that  the  terms  of  the  coefficient 
matrices  are  analytic  in  all  their  arguments.   It  remains  to 
be  shown  how  best  this  extension  can  be  made,  since  we  will 
only  be  interested  in  that  portion  of  the  solution  which 
lies  in  the  real  uv-plane. 

V7e  now  state  the  general  problem: 

We  wish  to  solve  the  complex  system  (1.16),  (1.17) 
where  X^  are  known  functions  of  u  and  v,  and  x  and  y  satisfy 
appropriate  initial  conditions.  We  allow  all  entities  to  take 
on  complex  values. 
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2.   TRANSONIC  CASCADES 

In  this  chapter,  we  use  physical  criteria  derived  from 
the  geometry  of  cascades  to  further  specify  the  mathematical 
problem,  especially  with  respect  to  the  form  of  the  singular 
solution  and  the  assignment  of  proper  initial  data.   This, 
in  turn,  leads  to  a  further  development  of  the  characteristic 
equations . 

1.    Cascade  Geometry 

A  cascade,  as  mentioned  before,  is  a  two  dimensional, 
periodic  array  of  blades.   Steady  flow  through  such  a 
structure  is  used  to  represent  the  regime  of  an  axial  flow 
compressor.   It  is  assumed  that  a  fluid  with  a  prescribed 
velocity  at  x  =  -<»  (called  the  inlet  velocity)  passes 
through  the  cascade  to  arrive  at  another  prescribed  velocity 

(called  the  exit  velocity)  at  x  =  +».   For  example,  in  Figure  5 

-+  .    -♦■ 

we  have  a  fluid  with  inlet  velocity  v,  and  exit  velocity  v,. 

-♦■ 
The  inlet  angle  a,  ,  is  the  angle  between  v^  and  the 

vertical    as  is  also  the  case  for   the  exit  angle  a-  and  v^- 

The  net  change  in  angle  (not  shown  in  the  figure)  undergone  by 

the  flow  from  x  =  -«=  to  x  =  +<»  is  called  the  turning  angle. 

The  gap  g  is  the  separation  between  blades  while  c  is 

the  chord  length.   In  practice  we  usually  use  the  gap-to-chord 

ratio  g/c,  or  its  inverse,  the  solidity  c/g,  which  measures 

the  relative  distance  between  blades.   We  can  think  of  the 


-19- 


thickness  t  as  being  the  maximum  blade  width  along  the  chord. 
Here  also,  we  are  usually  interested  in  the  relative  thick- 
ness, given  by  the  ratio  t/c . 

Last   (see  figure  )  ,   the  stagger  angle  B  specifies 
the  angle  at  which  the  blades  are  stacked.   For  our  purposes, 
vertical  stacking,  i.e.  6  =  90°,  will  almost  always  suffice. 
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2 .    Conjugate  Characteristic  Coordinates 

In  this  section,  we  use  the  theory  of  functions  of  a 
complex  variable,  with  special  reference  to  its  description 
of  incompressible  flow  about  the  unit  circle,  as  a  guide  in 
our  choice  of  characteristic  coordinates.   We  remark  that 
throughout,  we  will  be  concerned  with  flow  about  a  single 
body.   In  fact,  our  choice  of  coordinates  will  be  based  on 
compressible  flow  about  an  airfoil  ([  2  ] ,    pp.  15-19)  , 
the  periodicity  requirem.ents  of  cascade  geometry  being 
satisfied  later  as  a  result  of  the  form  of  the  singular 
solution. 

Clearly,  we  can  consider  all  quantities,  in  particular 
u  and  V,  to  be  functions  of  our  as  yet  undetermined  (complex) 
characteristic  coordinates  a  and  6.   It  is  in  the  choice  of 
these  coordinates  that  we  have  the  most  freedom,  since  any 
analytic  function  of  a  characteristic  coordinate  is  itself 
a  characteristic  coordinate.   We  will  exploit  this  freedom 
by  requiring  that  a  and  0  be  not  only  characteristic  coordi-- 
nates,  but  conjugate  characteristic  coordinates,  i.e.  that 

a  =  i 

whenever  u  and  v  are  both  real,  and  conversely.   This  will 
simplify  our  formulation  of  the  problem  and  reduce  the 
required  computation,  since  we  can  use  the  above  relationship 
to  identify  that  part  of  the  solution  which  lies  in  the  real 
uv-plane . 
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To  gain  an  accurate  understanding  of  the  role  of 
conjugate  characteristic  coordinates,  we  now  examine  the 
case  of  incompressible  flow  about  a  body.   More  specifi- 
cally, we  consider  the  two  dimensional,  steady,  irrotational 
flow  of  an  incompressible  inviscid  fluid  about  the  unit 
circle  in  the  xy-plane.   It  can  be  shown  ([24],  pp.  174-188), 
of  course,  that  the  complex  velocity 

—  i6 
(2.1)         w  =  qe    =  u(x,y)  -  iv(x,y) 

is  an  analytic  function  of  the  complex  variable 

z  =  X  +  iy  , 

where  u  and  v  are,  respectively,  the  x  and  y  components  of 
the  flow,  as  before.   This  result  follows  from  the  fact  that, 
in  this  illustration,  the  flow  equations  reduce  to  the  familiar 
Cauchy-Riemann  equations  for  u  and  -v, 

(2.2a)  u^  =  -Vy 

(2.2b)  u  =  V   . 

y     X 

Then,  if  we  assume  a  velocity  u  =  1  and  v  =  0  at  infinity, 
with  circulation  r  about  the  unit  circle,  we  can  describe  the 
flow  by  the  complex  potential 

(2.3)  x(z)  =  z  +  -  +  2ik  In  z  , 

where 

2tt 
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Moreover,  by  the  standard  theory,  we  know  that  x  is  analytic 
in  z,  and  further,  that 

dz        ^2     2 
Inverting  this  last  equation,  we  see  that 

(2.4)  z  = ^  ■■  I-) • 

(1-w-k  )-^^  +ik 

We  emphasize  here  that  the  relationship  between  z  and 
w  as  expressed  in  equation  (2.4)  has  a  dual  interpretation. 
The  familiar  one,  based  on  complex  function  theory,  is  that 
z  is  an  analytic  function  of  w,   and  hence  that  the  two 
variables  are  related  by  a  conformal  mapping.   For  reasons iO^f 
convenience,  we  introduce  the  com.plex  variable, 

(2.5)  C  =  (l-w-k^)-*-/^  +  ik  , 
so  that  we  have 

(2.6)  z  =  |- 
and, 

(2.7)  w  =  1  +  2ik5  -  ^    . 

In  effect,  then,  we  have  a  mapping  from  the  z  plane  to  the 
interior  of  the  unit  circle  in  the  C  plane,  and  then  to  the 

w  plane.   Figure  6  shows  such  a  mapping  for  k  =  1/2.  Observe 

2 
that,  as  a  result  of  the  branch  point  at  w  =  1-k  ,  the  mapping 

from  £;  to  w  is  double-valued.   Also,  z  has  a  simple  pole  at 

5=0  (or  w  =  1) .   We  remark  that,  up  to  this  point  in  our 
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illustration,  x,  y,  u  and  v  are  all  real  variables. 

On  the  other  hand,  we  may  introduce  the  hodograph 
transformation  and  carry  out  a  computation  similar  to  that 
of  Section  1.3  to  get  the  characteristic  equations 


(2.8a)  V  =  A  u 

a    +  a 


(2.8b)  x^  =  X^y^  , 

which  hold  on  the  characteristic  planes 

g(u,v)  =  constant, 
and 


(2.8c)  Vg  =  X_Ug 


(2.8d)  x^  =  X_Y^    , 


which  hold  on  the  characteristic  planes 

a(u,v)  =  constant. 


where 


X^  =  ±i 


Note  that  since  the  system  (2.2)  is  of  the  elliptic  type, 
the  characteristic  directions  X^  ,  and  characteristic  coordi- 
nates a  and  6,   must  be  complex.   In  addition,  x,  y,  u,  and  v 
are  now  assumed  to  be  complex  also,  as  a  result  of  analytic 
continuation  (see  Section  1.4).   Moreover,  it  can  be  shown 
that,  in  this  context,  w  and  the  complex  variable 


-24- 


(2.9)         w*  =  qe^^  =  u(x,y)  +  iv(x,y) 

a.ve  conjugate  characteristic  coordinates  ([2],  pp.  12-13,17),  i.e. 
we  may  set 

a  =  w 
* 

e  =  w 

in  ;2.P),   Cf  course  the  w  used  here  is  the  analytic  continua- 
tion of  that  introduced  in  equation  (2.1),  since  u  and  v  are 
now  complex.   Thus  we  may  also  view  equation  (2.4)  as  a 
relation  between  the  com.plex  variables  x  and  y,  and  the  conju- 
gate characteristic  coordinate  w.   We  note,  in  passing,  that 
for 


* 
z  =  X  -  iy 


v;e  have 


(1-w  -k  )  '   -  ik 
which  expresses  a  relaticfiship  between  x  and  y,  and  the  remaining 
coordinate  w  . 

Going  further,  we  may  also  carry  over  the  variable  E, 
as  defined  in  equation  (2.5),  subject  to  a  proper  analytic 
continuation,  and  in  a  similar  way  introduce 

(2.11)  n  =  (l-w*-k2)l/2  _  ^y.    ^ 

to  get 

*    1 

(2.12)  z      =   - 


n 


and 


(2.13)  w*  =  1  -  2ikn  -  n^ 
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Clearly,  E,   and  n  are  also  conjugate  characteristic  coordinates, 
since  they  are  analytic  functions  of  w  and  w  ,  and 

"*■       _ 

w  =  w  o  ^  =  n  . 

These  are  the  coordinates  that  we  use  in  practice,  since 
unlike  w  and  w  ,   their  relationships  to  z  and  z   are,  respec- 
tively, one-to-one  for  real  x  and  y. 

Let  us  assuiTie  therefore,  that  we  have  rewritten  the 
characteristic  equations  (2.8)  in  terms  of  E.   and  n .   We  can 
then  use  (2.7,  2.13)   to  solve  two  of  these  equations  (2.8a,c) 
explicitly,  since 

w+w* 


and 


^  =   2 


* 
w  -w 

V  = 


2i 


Thus  the  illustration  is  complete.   We  point  out,  however, 
that  although  in  this  case  we  also  know  x  and  y  in  terms 
of  C  and  n  as  a  result  of  (2.6),  (2.12),  this  is  not  generally 
true  in  the  transonic  case.   There,  we  will  have  to  solve  the 
equations  corresponding  to  (2.8b,d)  numerically. 

We  now  turn  to  compressible  flow  about  an  airfoil, 
carrying  our  formulation  over  in  an  analogous  way.   It  has 
been  shown  ([2],  pp.  21,22)  that  for  such  flows,  the  conjugate 
characteristic  coordinates  corresponding  to  w  and  w   are 

(2.14a)  s  =  h(q)e"^® 

and 
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(2.14b)  s*  =  h(q)e^®  , 


where 


(2.15)   h(q)  =  Kq[(c  -  /c^-c^  ) ^^^   ( /c^-q^  +  c) ] 


and  K  is  a  constant  such  that 


h(l)  =  1 


Continuing  our  analogy  to  the  incompressible  case,  we  generalize 
equations  (2.7),  (2.13)  and  define  C  and  n  by 


(2.16a)  s  =  C^^^ 


using  the  summation  notation,  and 


(2.16b)  s  =  C^n   , 


noting  that  the  complex  constants  C.  are  to  be  determined  later, 
It  is  clear  that,  since  s  and  s   are  conjugate  characteristic 
coordinates,  C  and  n /  as  defined  in  (2.16),  will  be  too. 

Thus  we  may  now  reformulate  the  characteristic  equations 
(1.16),  (1.17).   We  wish  to  solve 


(2.17a)  V  =  X_j^u 


(2.17b)  Xp  +  X_y   =  0 

on  the  curves  C   ,  given  by 

5(u,v)  =  constant, 
and 
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(2.17c)  V  =  X  u 

n   -  n 

(2.17d)  X   +  X  V   =  0  , 

on  the  curves  C_  ,  given  by 

ri(u,v)  =  constant  . 

We  remark  that  we  will  have,  in  effect,  solved  (2.17a,c) 
when  we  have  specified  the  C.  of  (2.16),  since  for  each 
E,   and  n  we  can  get  s  and  s   by  (2.16)  and  in  principle 
invert  (2.14)  to  find  q  and  G  and  hence  u  and  v.  Neverthe- 
less, we  must  still  solve  (2.17b,d)  numerically. 

Attempts  at  designing  transonic  airfoils  using  these 
equations  have  met  with  much  success.   It  must  be  remembered, 
however,  that  (2.16a,b)  were  simply  taken  over,  albeit  in  a 
more  general  form,  from  incompressible  flow  about  a  circular 
cylinder,  and  there  is  nothing  in  this  presentation  to 
guarantee  their  effectiveness.   In  fact,  it  was  found  neces- 
sary to  modify  them  to  obtain  the  results  of  this  paper. 
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3.    The  Singular  Solution 

In  this  section,  we  use  properties  of  the  flow  through 
a  cascade  to  work  out  the  form  of  the  singular  solution. 

Consider  the  complex  extensions  of  the  planar  coordi- 
nates X  and  y.   Our  imm.ediate  goal  is  to  write  them  in  the 
form 

X  =  Xg(c,n)  +  Xj^(?,n) 

y  =  Yg(5,n)  +  Yj^(5,n)  , 

where  X  and  Y„  are  regular,  and  X_  and  Y_  have  singularities 
appropriate  to  a  cascade. 

Toward  this  end,  we  first  let  (u  ,v  )  and  (^^''^a^ 
respectively  denote  the  inlet  and  exit  flow  velocities. 
Assuming  the  coefficients  C.  to  be  known,  we  can  determine 
complex  numbers  Cp,  and  E,      such  that 

-^®B        i 

(2.18a)         Sg  =  h(qg)e   "^  =   C^Cg 

and 

(2.18b)  s^  =  h(q^)e   ""  =  C^C^  , 


where. 


and 


%    =    ^^''l    '  «B  =  ^^^"'^^^' 


^l--l^<'  ^A  =  t-""'(^) 
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Supposing  the  flow  to  be  from  left  to  right,  and 
remembering  that  u  and  v  are  real  when  n  =  I ,  we  require 
a  solution  such  that: 

(1)  The  velocity  vector  (u,v)  -^    (u^,v  )  and  x  ^^  -<» 
in  the  real  physical  plane  as  (^,ri)  — >■  (Cgflg)« 

(2)  We  have  a  source  at  the  point  in  the  real  physical 
plane  corresponding  to  (^  ,^g). 

(3)  The  velocity  vector  (u,v)  -^   ^^a'^A^  ^^^   x  — >-  +<» 
in  the  real  physical  plane  as  {^,n)  — >-  (^a'^A^  ' 

(4)  We  have  a  sink  at  the  point  in  the  real  physical 
plane  corresponding  to  (S^f  ?,>)  • 

(5)  The  solution  is  multiple- valued,  i.e.  we  generate 
a  cascade  instead  of  a  single  blade. 

(6)  The  stream  function  ip  is  single-valued  for  any  real 
closed  path  around  any  single  blade. 

All  of  these  conditions  can  be  met  by  a  solution 
of  the  form 

(2.19a)   x(?,n)  =  x-'-(^,n)  In  (?-C^)  +  X^(C,n)  In  (C-?g)  +  x^(C,n) 

(2.19b)  y(C,n)  =  Y^(c,n)  in  ii-i^)   +  Y^(^,n)  in  (5-Cg)  +  Y^(c,n), 

where  X"""  and  Y"*",  i  =  1,2,3,  are  regular  functions   to  be 
determined  nuirierically .   Note  that,  though  we  have  been  able 
to  choose  characteristic  coordinates   such  that  u  and  v  are 
real  when  t\    =   \ ,    this  is  not  necessarily  true  for  x  and  y; 
their  dependence  on  ^  and  n /  unlike  that  of  u  and  v,  is  still 
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in  the  form  of  a  pair  of  unsolved  equations.   Hence,  after 
solving  these  equations   for  the  X  and  Y  ,  we  will  choose 

{2.20a)  y-^^^t    =   Re  [x{E,,l)  ] 

(2.20b)  y^.^^^^  =  Re  [y(C,a] 

as  our  real  solution. 
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4 .     The  Complex  Characteristic  Initial  Value  Problem 
for  Transonic  Flow  through  a  Cascade 

In  this  section,  we  use  the  expected  form  of  the  solu- 
tion to  derive  a  final  system  of  characteristic  equations, 
from  which  we  construct  three  complex  characteristic  initial 
value  problems  ,  one  for  each  of  the  three  pairs  of  functions 
X^  and  Y^,  i  =  1,2,3,  of  equations  (2.19).   This,  in  turn, 
leads  to  a  further  refinement  of  the  singular  solution. 

Consider  the  system  of  characteristic  equations  (2.17). 
We  recall  that  we  have,  in  essence,  already  solved  (2.17a,c), 
since  equations  (2  .14) -  (2  .  16)  have  the  effect  of  specifying 
u  and  V  as  functions  of  C  and  n .   It  follows  then,  that  we 
need  only  solve  the  remaining  equations.   Since  we  have 
already  determined  the  form  of  the  solution,  we 
substitute  (2.19)  into  (2.17b,d)  to  obtain 


x^  + 


A_y.  =  (X^  +  \_Y^)    In  (?-C^) 


X^+X   Y-"- 


2     2 

2     2  ^  "^^-^      3       3 

+  (x|+x_y|)  In  (C-?g)  +  -^r^ +  X^  +  X_Y^  =  0 

B 

on  n=  constant,  and 


+  X^  +  X^Y^^   =   0 
on  E,   =   constant. 
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We  can  satisfy  the  first  equation  by  setting 


(2.21a)  X^  +  X_yJ  =  0 

(2.21b)  X^  +  X_Y^  =  0 


and 

11  2     2 

7,                 r,         X^+X_Y^  X  +X  Y^ 

(2.21c)       X^  +  X  Y^  =  — +  ■= ^ 


on  n  =  constant.   Similarly,  we  can  satisfy  the  second  by 
requiring 

(2.21d)  X-"-  +  X  Y"'^  =  0 


(2.21e)  X^  +  X^Y^  =  0 


and 


n    +  n 

2      2 

r  +  X  Y^ 
n    +  n 


(2.21f)  x-^  +  X^Y^  =  0 

3     3 
on   C  =  constant.  Finally,  if  X  and  Y  are  to  be  regular 

solutions  of  (2.21c),  we  must  have 

(2.22a)  X-^  +  X_Y^  =  0 

on  C  =  5^  ,  and 

(2.22b)  X^  +  X  Y^  =  0 


on  5  =  5  .   Clearly  we  can  use  the  equations  (2.21)  to  formu- 
late  an  initial  value  problem  for  each  of  the  three  pairs  of 
functions  X  ,Y  ,  i  =  1,2,3  once  we  have  specified  appropriate 
initial  data  ([2l]f  PP.  9-12). 
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To  do  that,  we  consider  the  two  functions  X  and  Y 
and  their  corresponding  equations  (2.21a,d).   For  a  properly 
posed  characteristic  initial  value  problem,  we  must  prescribe 
initial  data  on  one  characteristic  from  each  family.   We  can 
construct  the  data  to  within  an  additive  constant  on  the 
characteristic  plane  E,   =    E,        by  integrating  (2.21d)   in 
that  plane,  subject  to  the  condition  (2.22a).   Thus  we  have 
only  to  pose  initial  data  on  some  n  characteristic,  say 
n  =  n,-   We  will  use  the  value  n,  =  E,      in  hope  of  exploiting 
symmetry.   Ideally,  since  we  are  only  interested  in  those 
points  for  which  n  =  C  ,  we  would  like  the  initial  data  on 
n  =  n,  to  satisfy 


x^(5,n^)  =  x^(5,i^)  =  x^(C^,!) 


and 


Y^(C,n^)  =  Y^(?,?^)  =  Y^(?^,i) 


We  could  then  use  the  Schwarz  Reflection  Principle  to  signi- 
ficantly reduce  our  computation.   Unfortunately,  this  requires 
that  X  (^  ,f-)  and  Y  (E,    ,1,    )      be  real,  which,  in  general, 

they  are  not.   To  overcome  this  difficulty,  we  define  two 

^1  ^1 

new  functions  X  (C/H)  and  Y  (Cn)  by 

x^(C,n)  =  x-^(?,n)  +  ix^ 

Y^(5,n)  =  Y^(C,n)  +  iY^  , 

where  X,  and  Y,  are  real  constants  to  be  determined  later, 
A      A 

^1      ^1 
and  X  and  Y  have  the  desired  symmetry  properties  m  their 

initial  data. 
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A  similar  procedure  can  be  followed  in  posing  the  initial 

2       2 
value  problem  for  X   and  Y   as  solutions  of  (2.21b,e)  with 


B 


initial  data  prescribed  on  the  characteristic  planes  C  =  C 

and  n  =  Ht,  =  Co   in  terms  of  the  symmetric  functions  X  (^,ri) 
B     B 

/^2 

and  Y  {^,n)  and  the  real  constants  Xg  and  Y^. 

3        1 
For  X  and  Y  ,  we  choose  yet  another  pair  of  character- 
istics, say  C  =  5   and  n  =  ip  =  I^  which  are  far  enough  away 
from  E,.    and  ?„  ,  and  n,  and  n„  ,  respectively,  to  prevent 

Ac  A        rJ 


computational  error.  The  choice  of  initial  data  for  this 

'c'-c 


problem  will  be  discussed  later.   We  can  assume  X  (5_,^p) 


and  Y    (E,    ,1    )      to  be  real,  since  the  expected  solution  uses 
only  their  real  parts,  and  therefore  assign  symmetric 
initial  data  directly. 

Thus  we  have  a  characteristic  initial  value  problem 
in  three  parts : 
(I)    Find  X-^{^,ri)  and  Y^{E,,t\)    satisfying 

(2.23a)  X^  +  X_Y^  =  0 

on  n  =  constant,  and 

(2.23b)  9}   +  X_^Y-'-  =  0 

on  C  =  constant,  with  initial  data  X  (Sj.,n)  and  Y  (C,»n) 
obtained,  up  to  an  additive  real  constant,  by   solving 
(2.2  3b)  on  C  =  Cjv  f  subject  to  the  condition 

(2.23c)         (X-'-+iX^)  +  X_(Y-'-+iY^)  =  0  . 
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The  initial  data  O"^  n  =  n^  =  ^a   ^r®  given  by 


(2.23d)  X-^(5,n^)  =  X^(5^,!)  , 

and 


(2.23e)  Y^(C,n^)  =  Y^(C^,U  . 

(II)   Find  X^(5,n)  and  Y^(C,n)  satisfying 
(2.23f)  X^  +  A_Y^  =  0 

on  ri  =  constant,  and 

(2.23g)  X^  +  A^Y^  =  0 

^2  ^2 

on  5  =  constant,  with  initial  data  X  (?^,n)  and  Y  {^^,r\) 

obtained,  up  to  an  additive  real  constant,  by  solving 
(2.23g)  on  5  =  5d  '  subject  to  the  condition 

(2.23h)         (X^+iX„)  +  X  (Y^+iY„)  =  0  . 

B     ~      a 

The  initial  data  on  n  =  n^  =  Cr.  ^^^  given  by 


(2.23i)  X^(5,i1g)  =  X^(Cg,f)  , 

and 


(2.23J)  Y^{K,r]^)    =  Y^(Cg,!)  . 


(Ill)  Find   X^(C,n)  and  Y'^(C,n)  satisfying 

112     2 
^       ^    X-^+A_Y     X  +A  Y 

(2.23k)       X^  +  X  Y^  =  -p T^ +  -r f 


on   n  =  constant,  and 
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(2.23i)  X-^  +  \_^Y^   =  0 

3  3 

on  C  =  constant,  with  initial  data  X  (C_,ri)  and  Y  {E,    ,r]) 

obtained,  up  to  cin  additive  real  constant,  by  solving 
(2.23£)  for  Y^   on   E,   =   i^   where 

(2.23m)  y?{E,^,T\)    =  G(n)  , 

and  G  is  an  as  yet  undetermined  analytic  function.   For  the 
initial  data  on  n  =  Hp  =  fp  /  we  set 


(2.23n)  y?(E,,T\^)    =  Y?{E,^,l)     , 

and 


(2.230)  Y^(C,n^)  =  Y^iK^.l) 
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5.     The  Arbitrary  Constants  in  the  Singular  Solution 

In  this  section,  we  use  previously  specified  properties 
of  the  solution  to  determine  the  unknown  constants  X   ,  X^  , 
Yj.  and  Y^  ,   and,  in  the  process,  we  explain  the  role  they 
play  in  the  geometry  of  a  given  cascade. 

Consider  a  simple  closed  curve  C  having  a  counter- 
clockwise orientation,  which  lies  inside  the  body,  contains 
the  two  points  (C,,Ca)  arid  {^^,E^)    in  its  interior,  and  is 

A   A  si       a 

Otherwise  arbitrary  (see  Figure  7  ) .   The  image  of  C 
in  complex  xy-space  under  the  mapping  given  by  the  solution 
of  our  problem  will  be  a  closed  curve  in  the  real  physical 
plane,  say  F,  which  contains  a  single  blade  in  its  interior. 
Note  that  the  exterior  of  C  will  map  onto  the  interior  of  T, 
and  that  T   will  have  clockwise  orientation. 

Substituting  X'^'+iX^  ,  X^+iXg  ,  Y'^'+iY^  ,  and  Y^+iY^  for 

12   1       2 
X  ,X  ,Y  and  Y   in  (2.19),  respectively,  we  can  show  that 


(2.24a) 
and 


dx  =  2Tr(X^  +  Xg)  , 


(2.24b) 
which  implies  that 


({)  dy  =  2Tr(Y^  +  Y^)  , 


(2.25a) 

and 

(2.25b) 


^A  ^  ^B  =  0 


^A  ^  ^B  =  °  ' 
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since  we  require  the  solution  to  be  single-valued  about 
such  a  curve. 

Similarly,  by  requiring  the  stream  function  ^   to  be 
single-valued  on  r,  we  get 


(2.26)   A  dil-  =  j)  pu  dy  -  pv  dx  =  27T  [p^u^Y^+  PgU^Yg-  P^v^X^ 

r     r 

which,  using  (2,25),  implies  that 


(2.27)         [Pa"a-Pb^B^^A  =  (Pa^A-Pb^B^^A  ' 

where  we  assume  that  p^^  and  p^   are  already  known  in  terms 
of  the  velocity  components. 

Observe  that,  had  we  chosen  curves   C,  or  C_ 
(see  figure)  containing  (C.'^a^  alone  or  (^  ,|  )  alone  in 
the  interior,  respectively,  we  would  have  obtained  the  curves 
r^  or  T^    such  that 


(2.28a)  (p  dx  =  2tt  X^  , 

^1 


(2.28b)  i)  dy  =  2Tr  Y^  , 

^1 

(2.28c)  (D  dx  =  2tt  Xg  =  -2tt  X^  , 

and  2 


(2.28d)  d)  dy  =  2tt  Yg  =  -In   Y^  . 

^2 

Thus  we  see  that  integration  along  either  of  these 

curves  increases   the  x  coordinate  by  +  X  and  the  y  coordi- 
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nate  by  +  Y.  which,  in  the  physical  plane,  corresponds  to 
jiomping  to  another  blade  of  the  cascade  (see  Figure  5  )  . 

Based  on  the  above  discussion,  we  define  the  stagger 
angle  3  and  gap  g  by 

-1  ^A 
(2.29a)  3  =  tan  -^  TT^  , 

^A 

and 

(2.29b)  g^=  X^  +  Y^   . 

Recall  that  for  g  to  be  a  useful  description  of  the 
blade  separation,  it  must  be  normalized  in  such  a  way  that 
it  measures  the  distance  between  blades  relative  to  some 
particular  blade  parameter,  say  the  chord  length.   Unfortunate- 
ly the  blade  itself  an^  consequently  its  chord  length  is 
not  known  until  after  the  computation  is  finished.   Hence  we 
use  instead  the  relation 

(2.30)  g  =  Ua-^b'""''  • 

The  reason  for  this  somewhat  imperfect  choice  is  that, 
other  parameters  remaining  constant  as  E,      ^   E,      ,    our  entire 
model  will  approach  the  airfoil  case.   Equation  (2.30)  is 
consistent  with  that  occurrence  in  that  it  says  the  gap 
becomes  infinite,  i.e.  we  have  a  single  blade. 

If,  in  addition,  we  require  B  =  90°,  which  means  that  the 
blades  are  stacked  vertically,  we  can  use  (2.29)  and  (2.30)  to 
find  X  and  Y   ,  and,  after  prescribing  three  of  the  four  inlet 
and  exit  velocity  parameters,  use  (2.25)  and  (2.27)  to  deter- 
mine all  the  remaining  constants. 
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6.    The  Initial  Data 

In  this  section, we  present  the  general  form  of  the 
initial  function  and  describe  the  procedure  by  which  we 
determine  its  unknown  constants. 

Recall  that  the  initial  data  for  esch  of  the 
functions   X  ,  Y  ,  X   and  Y   have  already  been  determined 
to  within  an  additive  real  constant  in  the  previous  two 
sections.   It  remains  for  us  to  prescribe  an  analytic 
function  G(n)  to  be  used  as  initial  data  for  X  (^p^n) 

in  the  condition  {2.23m).  Y  (^p/l)  will  then  be  found  by 

-      3  3 

integrating  (2.23£)  and,  since  n^  =  Cp  »  X  {K,r\^)    and  Y  (C,n^) 

will  be  determined  by  the  symmetry  conditions (2 . 23n,o) . 

Once  we  have  determined  G,  we  compute  the  solution  and 
plot  the  streamline  4*  =  0  in  the  physical  plane.   When  ^   has 
been  properly  initialized,  this  curve  will  be  the  two  dimen- 
sional cross  section  of  a  single  turbine  blade.   Hence,  we 
will  refer  to  the  streamline  4;  =  0  in  the  physical  plane, 
as  well  as  in  the  C  =  ii  plane,  as  the  body.   All  other  blades 
in  the  cascade  will,  of  course,  have  exactly  the  same  shape. 
Observe  that,  as  we  have   said  before,  the  body  will  not,  in 
general,  look  like  a  turbine  blade.  The  basic  plan  then,  is 
to  pick  G(n)  so  that  it  does. 

Theoretically,  we  have  a  considerable  amount  of  freedom 
in  our  choice  of  initial   function,  since  the  only  restriction 
on  G  is  that  it  be  analytic.   In  practice,  however,  we  proceed  in 
a  way  similar  to  that  used  in  the  airfoil  case  ([2],  p.  19)  and  set 
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M  N  . 

(2.31)   G(n)  =  I    (h-^+    iA   )  In  (n-B.  -iB   )+  I       {h.^+ih.^)^^, 
j  =  l   -'      -'  -*     -'    j=M+l   -'     -^ 

where  the  A.„  and  B.„  ,  i   =  1,2,   are  arbitrary  constants. 
G  will  be  completely  determined  when  we  have  prescribed  these 
2N+2M  real  constants . 

Note  that  once  we  have  chosen  the  B.p  ,  G  will  have  a 
logarithmic  singularity  at  each  of  the  points 

n  =  B.^  +  iB.2  , 

which  of  course  lie   outside  the  region  of  flow. 
These  singularities,  which  play  the  role  of  sources  and 
sinks  in  the  initial  plane  5  =  C^  ,  are  used,  assuming  the 
A.„  to  be  small,  to  make  local  adjustments  to  the  streamline 
representing  the  body.   Our  procedure  will  be  to  choose  the 
B.»  and  then  use  physical  constraints  and  automation  paths, 
i.e.  polygonal  curves  in  the  ?  =  n  plane  along  which  we  expect 
the  body  to  pass,  to  determine  the  remaining  constants  in  G. 
Both  the  B . p  and  the  automation  paths  are  chosen  by  a  combina- 
tion  of  past  experience,  and  trial  and  error. 

Our  physical  constraints  fall  into  two  categories: 
those  which  apply  to  the  nose,  and  those  which  apply  to  the 
tail  ( [  2],  pp.  59-61) . 
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(A)    The  Nose  Constraints 

Each  blade  will  have  a  stagnation  point  near  the  leading 
edge  which  we  call  its  nose  (see  Figure  5  ) .   The  image  of  that 
point  in  the   £;  =  n   plane  will  have  coordinates  (Cjjffjj)  / 
where  E,      satisfies 

(2.32)  C^Cj  =  0  . 

This  follows  from  equations  (2.14,16)  and  the  fact  that  h 
is  zero  when  q  is  zero.   Thus  we  see  that  we  can  determine  i 
once  the  C  are  known,  though  we  may  still  have  to  choose 
between  several  roots. 

Since  the  body  passes  through  the  nose,  we  set 

where  the  superscript  N  indicates  that  i'    is  to  be  evaluated 

at  (5  ,C^j)  .   Initializing  ^   in  this  way  has  the  effect  of 

choosing  the  streamline  iJj  =  0  as  the  body.   The  value  of  ii   at 

any  other  point  will  be  found  by  integration. 

The  nose  constraints  arise  from  the  fact  that  the  values 

of  X  and  x  at  the  nose  depend  on  R^  ,  the  radius  of  curvature 
u      V  ^        C 

of  the  blade  surface  at  the  nose,  and  a,  the  angle  between 
the   inlet  velocity  vector  and  the  direction  of  the  central 
streamline  at  the  nose.   Indeed,  it  can  be  shown  that 

3x  N     1 
(2.33a)  (|^)   =  K-^(Rj^,a) 

(2.33b)  (|^)   =  K^(R^,a)  , 
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1       2 
where  K  and  K  are  known  functions.    Clearly,  then, 

these  equations,  which  will  be  our  nose  constraints,  can 

be  introduced  as  a  matter  of  convenience,  ir.aking  it 

possible  to  specify  R^  and  a  as  inputs. 

(B)    The  Tail  Constraints 

Unlike  that  of  the  nose,  (C_,Cm)/  the  image  of  the 
tail  in  the  ^  =  n  plane,  is  not  determined  by  equations 
(2.14),  since  the  flow  velocity  at  that  point  is  not  zero. 
Instead,  E,^   is  specified  as  an  input  parameter  of  the 
program,  and  its  final  value  is  determined  by  trial  and  error. 
We  rem.ark  that,  assum.ing  we  have  chosen  a  finite  value  of  E, 
which  does  in  fact  locate  the  tail,  equations  (2.14)  imply 
that  the  tail  velocity  is  finite  also. 

Our  first  constraint  is  that  the  body  go  through  the 

tail,  that  is,  that 

T 
ij^   =  0  , 

where  the  superscript  T  indicates  that  ^   is  to  be  evaluated 
at  (5^,1^). 

Second,  sc  tha<-  vre  in  f&ct  have  a  tail,  vve  require 

Since 

32  _  az  9w 
8C    3w  3C  ' 
and 

3_z  _  3x      3y 
3w  -  8u      8u  ' 
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we  see  that  we  will  have  achieved  this  result  if 

T       T 

(|^)   =  (|^)   =  0  . 
3u      9u 

Thus  we  have  the  three  tail  constraints 

(2.34a)  ip"^  =  0 

-   T 
(2.34b)  (|i)   =  0 

(2.34c)  (|^)   =  0  . 

Conditions  (2.34b,c)  constitute  our  formulation  of  the 
Kutta  condition. 


(C)    The  Automation  Paths 

Finally,  we  can  prescribe  certain  paths  in  the  subsonic 
region  of  the   C  =  n   plane  which  we  call  automation  paths, 
through  which  we  want  the  body  to  pass  ([2  ]/  PP-  36-39). 
Choosing  a  mesh  along  each  pach,  we  obtain  the  points  (Cu/Ci,)  / 
k  =  l,r,    at  which  we  desire 

Although  we  will  see  below  that  it  is  usually  impossible 
to  satisfy  this  condition   for  every  index  k,  we  can  still 
use  it  to  justify  a  least  square   fit.   In  substance,  choosing 
automation  paths  is  a  more  convenient  way  of  prescribing 
initial  data  than  directly  picking  the  remaining  constants  in  G. 
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We  now  use  the  above  three  categories  of  constraints 
to  determine  the  A.,  in  (2.31).   We  assume  that  the  B.„ 
have  already  been  chosen. 

Let   A, ,  =  1  and  set  all  the  other  A.„  =0.   Using 

11  jl  ^ 

this  initial  function,  we  determine  a  solution  and  its 

N  NT  T 

corresponding  values  (9x,/9v)  ,  (3x,/9u)  ,  ii,       ,     (9x^/3u)  , 

T       k 
(9y,/9u)   and  ^,       ,  k  =  l,r.   Repeating  this  procedure  for 

N  NT 

each  A..  ,  we  obtain  the  values  (9x./9v)  ,  (Sx./9u)  ,  tp.   , 

(9x./9u)^,  (9y./9u)'^  and  ^  ^^    ,    k  =  l,r,   where   i  =  l,n 

and  n  =  2N.   But,  by  the  principle  of  superposition,  any 

linear  combination  of  solutions  will  again  be  a  solution. 

In  fact,  we  assume  that  we  have  constructed  such  a  solution 

using  the  as  yet  unknown  constants  c  ,   i.e.  for  the 

solutions  xi  ,  y.  ,  i  =  l,n  we  have 

n 
X  =   y   c .  X. 

i=l   ^  ^ 

and 

Our  nose  and  tail  constraints  then  become: 

n      9x.  N     , 
(2.35a)  I      c.   (^)   =  K-^(R_,a) 

.  ''  T    1  oU  C 

1  =  1 

n      9x.  „     -J 
(2.35b)  y   c.  (^)'^  =  K^(R^,a) 

•  ^-1    1  OV  I- 

r        T 
(2.35c)  I       C.\lj.^  =  0 

i=l   ^  ^ 
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n     3x.  T 
(2.35d)  I      c.(g^)   =  0 


i=l 


n     8y .  T 
{2.35e)  I      c.  {^)   =  0  , 

i=l   ^  '^^ 


to  which  we  add  the  normalization, 


n 
(2.35f)  y   c.   =   1 

i=l   ^ 


We   now  introduce  a  change  of  interpretation,  viz. 
we  view  (2.35)  as  a  system  of  equations  to  be  solved  for 
the  c..   Obtaining  them  in  this  manner  is  equivalent  to 
satisfying  the  constraints  by  a  proper  choice  of  the  A . „ . 

At  this  point,  of  course,  since  we  have  six  equations 
in  n  unknowns,  where  six  is,  in  general,  less  than  n,  our 
system  is  underdetermined.  Ideally  we  would  like  to  remedy 
the  situation  by  adding  the  automation  conditions 

r     k 
(2.36)  I      ci).    =  0  ,      k  =  l,r  . 

i=l   ^  ^ 

Unfortunately,  this  would  leave  us  with  6+r  equations  in 
n  unknowns,  and,  since  6+r  is  usually  greater  than  n,  the 
system  would  be  overdeter mined.   For  this  reason,  we  do 
not   use  (2.36),  requiring  instead  that  the  body  fit  the 
automation  paths  in  a  least  square   sense,  subject  to  the 
constraints  expressed  in  (2.35). 

Our  final  statement  of  this  problem  will  be  in  matrix 
form.   Let  the  equations  (2.35)  be  written 
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BC  =  F  , 

where  C  is  the  cclumn  vector  containing  the  c  ,  and  B  and 
F  have  been  chosen  accordingly.   We  then  define 


'^■o  =  {Jr'^'} 


where  the  oj  •  are  suitably  chosen  weights,  and  we  pick  the 
matrix  A  so  that  we  may  write  equation  (2.36)  in  the  form 


AC  =  0  . 
Thus,  we  may  restate  the  problem:   Minimize  the  expression 

subject  to  the  constraint 

BC  =  F  . 

Note  that  lACl   refers  to  the  Euclidean  norm  of  that  matrix, 
and  that  x   is  an  input  parameter  which  is   used  to  weight 
the  importance  of  the  least  square   fit  relative  to  minimizing 
the  size  of  the  c.   For  example,   x  =  0   tells  the  program 
to  minimize  the  c.  without  attempting  any  fit  at  all. 

Using  Lagrange  multipliers,  the  problem  reduces  to 
solving  the  system 

BC  =  F  , 

(x^a'^A  +  Iq)C  +  B^X   =   0  , 
where 
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(Jj 


U), 


0 
2 


I„  = 


0   (i) 


and   X   is  the  column  vector  of  multipliers.   Here,  of 

XT  .         . 

course,  A  and  B   denote  the  respective  matrix  transposes, 
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3.   THE  COMPUTATION:   A  CONCEPTUAL  OVERVIEW 

In  this  chapter  we  give  a  descriptive  account  of  the 
numerical  method  of  solution,  emphasizing  the  basic  strategy 
of  the  computer  program,  as  opposed  to  the  details  of  its  use. 

1.    The  Basic  Diagram 

As  a  starting  point  for  the  discussion,  we  first  take  a 
preliminary  look  at  a  typical  example  of  the  graphic  output 
of  the  program  (see  Figure  8  ) . 

Figure  8A  contains  the  computed  blade  profile,  indicating 
the  characteristics  (or  Mach  lines)  of  its  supersonic  zone,  and 
showing  above  the  Mach  niimber  distribution  along  the  blade 
surface. 

Figure   8B  ,   though  plotted  on  a  single  pair  of 
coordinate  axes,  variously  represents  either  the  initial  plane 
n  =  n   ,  or  the  n  =  ^  plane*  which  we  shall  loosely  refer  to 
as  the  hodograph  plane.   This  dual  representation  is 
accomplished  in  the  following  way:   the  position  of  a  given 
point  relative  to   the  coordinate  axes  of  the  diagram  deter- 
mines its  E,   coordinate.   Its  n  coordinate  is  then  given  by 
n  =  n   or  n  =  C  /  depending  on  which  plane  is  under  discussion. 
The  reader  is  forewarned  that  while  some  elements  of  the  graph 
belong  to  the  initial  plane  and  some  to  the  hodograph  plane, 
there  are  elements  that  belong  to  both  planes,  depending  on 
the  interpretation.  In  fact,  subject  to  proper  interpretation. 
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this  figure  depicts  any  plane  in  which  n  =  constant. 

The  basic  Figure   8B   is  composed  of  the  following 
elements: 

(i)    The  sonic  line/locus 

The  curve  C  is  defined  to  be  the  set  of  points  at  which 

(3.1)  iC^^^I  =  |h(q)e"^®|  =  h(c^)  . 

Note   that   in   the  hodograph  plane , 

_  * 

C   =   ri°*s   =   s      =»  h(q)       is    real. 

Thus,  equations  (2 . 14 ) -  (2 . 16)  and  (3.1)  imply  that 

q  =  c^  , 

i.e.  the  flow  is  exactly  sonic;  hence  we  call  C  the  sonic 
line.   Observe  that  though  the  sonic  line  contains  the  image 
of  the  real  subsonic  regime  in  its  interior,  the  supersonic 
regime  maps  onto  another  surface  entirely  (see  Section  3.4). 
Suppose  now  we  consider  C  to  be  a  curve  in  the  initial 
plane  n  =  n   /  whose  points,  however,  still  satisfy  (3.1). 
Here,  of  course,  h(q)  is  complex.   Nevertheless,  the 
E,   =  constant  characteristic  through  each  point  of  C 
terminates  on  the  sonic  line,  i.e.  the  points  of  C  in  the 
n  =  n   plane  are  mapped  into  the  sonic  line  by  the  character- 
istic curves.   For  this  reason,  we  call  it  the  upper  sonic  locus, 
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the  analogous  curve  on  C  =   ?,  t)6ing  called  the  lower  sonic 
locus. 

To  Slim  up,  the  curve  C  is  called  the  sonic  line  or 
the  sonic  locus,  depending  on  whether  we  choose  n  =  f 

or  n  =  n^,. 

(ii)   The  body 

This  pattern  of  points,  indicated  in  Figure  8B, 
represents  the  ij;  =  0  streamline  in  the  hodograph  plane. 
We  remark  in  passing,  that  the  subsonic  region  of  flow, 
which  is  obviously  outside  the  blade  in  the  physical  plane, 
maps  onto  the  interior  of  the  body  image  shown  here. 

(iii)  Stagnation  points  and  branch  points 

The  stagnation  points  are  the  roots  of 

(3.2)  C^C^  =  0. 

The  branch  points  are  the  roots  of 

(3.3)  1^  (C^C"-)  =  0  . 

Their  significance  will  be  discussed  later. 

(iv)   The  logarithmic  singularities  of  the  initial  function 

In  Section  6,  we  saw  that  the  initial  function  G(ri) 
has  logarithmic  singularities  at 

-52- 


(see  equation  (2.31)),  or  equivalently,  that  these  points 

are  the  locations  of  the  singularities  of  X  and  Y  in  the 

.  .  .  3      3 

initial  plane  E,   =   E.    .      Also,  by  reflection,  X  and  Y  have 

similar  singularities  at 


in  the  n  =  n   =  f  plane.   In  practice,  we  input  the 
^     c 

latter  set  of  singularities,  which  are  shown  in  the  figure, 
and  determine  the  former  by  reflection. 

There  are  two  basic  types  of  integration  paths: 

(v)    The  subsonic  paths 

These  paths  are  entirely  contained  inside  the  sonic 

line/locus  and,  aside  from  those  used  to  search  for  the  body, 
include 

(a)  the  constraint  paths  -  two  paths  which  go  to  the  nose 

and  tail  and  are  used  to  compute 
the  respective  constraint  condi- 
tions (see  Section  2.6) 

(b)  the  automation  paths  -  paths  which  indicate  the  expected 

location  of  the  body  in  the  hodo- 
graph  plane  (see  Section  2.6). 

(vi)   The  supersonic  paths 

These  paths  go  from  the  initial  ^-characteristic  to  the 
image  of  the  real  supersonic  regime  in  (5,ri)  space  (see  Section 
3.4)  . 
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2.     Subsonic  Integration 

For  the  moment,  we  turn  out  attention  to  the  simple 
polygonal  path  outlined  in  dark  ink  in  the  subsonic  region 
of  Figure  8C.    We  will  see  how  this  curve  is  used  to 
construct  a  path  in  the  hodograph  plane,  along  which  we 
search  for  points  of  the  body,  or   equivalently ,  zeros 
of  the  stream  function. 

As  a  first  step,  we  consider  this  path  to  be  in  the 
initial  plane  n  =  n  •   It  consists  of  line  segments  connecting 
the  points  E,^    ,    E.^    ,    E.^    ,    ^^    ,    E,^   and  Cp  in  that  order.   From 
this  information,  we  can  construct  its  conjugate  path,  viz. 

the  path  connecting  "Ha  '  '^b  '  ^C  '  ^D  '  ^E  '  ^'^^  ^F  '  ^^^^^ 
each  n  is  the  complex  conjugate  of  its  corresponding  E, ,    and 
where  we  assume  the  path  to  be  in  the  E,   =   E,      =  n   initial 
plane.   Now,  we  choose  a  net  of  points  along  each  path. 


and 


where 


^A  ~   ^l'^2""  '^M        ^F  ' 


"^A  "  r^^,T]^,...,r^^   =  rip 


n^  =  Cj^  ,        i  =  1,2,  .  .  .  ,M. 


If  we  consider  all  possible  pairs  of  points  in  this  mesh, 

U^^r].)     ,      i,  j  =  1,2,.  .  .  ,M, 

we  obtain  the  two-dimensional  array  of  complex  numbers 
outlined  schematically  in  Figure  9. 
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It  is  on  this  grid  that  we  solve  (see  next  section  for  finite 
difference  procedures)  our  threefold  initial  value  problem 
numerically. 

Clearly,  we  can  use  equations  (2.23a-e)  to  generate 

.  .  ■  ^1     ^1 

initial  data  for  X  and  Y  on  C  =  C,  and  n  =  n,  .   Thus  we 

know  those  functions  up  to  an  additive  real  constant  in  the 
first  column  and  last  row  of  the  array.   We  then  use  equa- 
tions (2.2  3a,b)  to  find  their  values  at  all  the  remaining 
points.   A  similar  procedure  works  for  X  and  Y   in  their 
region  of  definition,  where  we  generate  initial  data  on 

C  =  ?g  and  n  =  rig. 

3       3 
We  now  consider  X   and  Y   m  the  quadrant  defined  by 

their  initial  characteristics  C  =  C   and  n  =  n  .  Again,  we 

know  these  functions  in  their  initial  planes,  having  already 

determined  G(ri)  (see  equations  (2.23k-o)  and  Section  2.6). 

"^  i     ^  i 
Since  we  already  know  X   and  Y  ,  i  =  1,2,   everywhere,  it 

112      2 
follows  that  we  also  know  X  ,  Y  ,  X  and  Y  everywhere,  and 

can  solve  the  inhomogeneous  equations  (2.23k,l)  for  X  and  Y 

over  the  entire  grid.   Thus  we  see  that  we  can  determine  the 

unknown  functions  at  every  point  of  the  grid  for  which  they 

are  defined.   Once  this  is  done,  all  other  functions,  e.g. 

i|;,p,  etc.   can  be  determined  also. 

Consider  now  the  diagonal  entries  of  the  array.  Since 

the  paths  are  conjugate,  we  have 

U^,r]^)    =    (Ci,q)  ,        i  =  1,2, ...,M, 
i.e.  each  point  on  the  diagonal  is  in  the  real  hodograph  plane, 
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Hence,  for  each  pair  of  conjugate  initial  paths,  we  obtain 
a  solution  along  a  single  path  in  the  real  hodograph.  Any 
stream  function  zeros  on  the  diagonal  will  map  into  the  body 
in  the  real  physical  plane.   Of  course,  choosing  good  integra- 
tion paths  (e.g.  paths  which  intersect  the  body)  is  the  task 
of  the  user,  and,  in  the  case  of  automation  paths,  requires 
some  skill. 

In  addition,  we  point  out  that  we  can  eliminate  all  our 
computation  below  the  diagonal,  since  e.g. 


x-'-(q,nj)  =  x^(q,nj)  =  x^(n^,Cj)  . 

Similar  symmetry  conditions  for  the  other  functions  allow  us 
to  reduce  the  amount  of  calculation  by  almost  a  half. 

We  remark  that,  except  in  the  case  of  the  constraint 
and  automation  paths,  we  usually  prescribe  integration  paths 
as  a  fork  with  many  stems  (see  Figure  8B ) .    This  allows  us 
to  pick  up  more  than  one  body  point  in  a  region  without  having 
to  repeat  the  entire  integration. 
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3.    The  Finite  Difference  Scheme 

Consider  the  grid  of  an  arbitrary  path  (see  Figure  10) 
in  which  we  have  used  the  shorthand  notation 


P.  .  =  (q,nj)  . 

We  assume  that  values  for  X  and  Y  at  P^ ,  have  been 
prescribed,  and  will  see  how  the  program  determines  these 

functions  at  other  points  of  the  grid. 

^l  ^1 

First,  we  need  to  know  X  and  Y  on  their  initial 

C-plane,  where  n  =  n  .  ,  j  =  1,2,...,M,  i.e.   we  wisb  to  solve 

equation  (2.2  3b),  subject  to  the  condition  (2.23c),  where 

all  functions  are  to  be  evaluated  on   E,   =   i      (see  Section  2.4). 

To  do  this,  we  proceed  point  by  point  up  the  first  coliomn 

of  the  grid. 

As  an  example,  consider  the  point  Pi2-   ^^®  already  know 

the  solution  at  P, ,  ,  and  assiune  the  approximations 

^1  ^  X^(Pi2)-^'(Pll) 


n      ^2  "  "^1 


and  ^,      ., 

Y  (P   )-Y  (P   ) 
-1  ^  ^  ^^12^  ^  ^^11^ 


n        ^2  "  "^1 


Substituting  in  (2.23b),  we  get 


X^  (Pl2^  "^^  ^^11^      *  f^^  ^^12^  "^'''  ^^11^ 


+  X  . 

^2-^1  "- 


^2  -  ^1 


=  0  , 


which  becomes  our  first  difference  equation 
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(3.4a)   X^(P^2)-^'''(Pii)  +  ^*(^'''(Pl2^~^'''^^ll^)  =  °  • 


Here,  we  have  used 


X.  = 


X^(P,l)-fX^(P^2^ 


where  we  determine  both   X^(P,2)  iri  the  following  way. 
Clearly, 


and 


implying  that 


h(q(P^2)  )e 
h(q(P^2)^^ 


-ie(P   ) 

=  c  f: 


^i^2  ' 


h^(q(P^2^)  =  t"^i^i^  f^i'^2^  ' 


which  we  invert  numerically  to  find  qlP^^).   Similarly,  we 
obtain  e(P^2)  ^Y  solving 

216  (P^2'    5^r,^ 


^i^l 


to  get 


(P12)  =  21  1" 


c.n^ 


c,,l 


Knowing   q(P^_)   and   e(P^„),  we  can  find  cCP,-)*  u(P,_), 
cind  v(P,2)  t>y  way  of  the  standard  equations  of  fluid  dynamics, 
e.g.  Bernoulli's  law,  etc.   We  then  determine   X^(P  „)  ^Y 
substitution  into  equation  (1.15). 

Vie   c,et  cur  second  equation  by  evaluating  (2.23c)  at  P12' 
which  yields 
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(3.4b)    X^(P^2)  "^  ^^A  "^  A_(Pi2^  '^^^^12^  ^    ^^p}       "   ° 


where  X.,  ^_(Pi2)  ^'^d  ^   are  already  known. 

We  can  now  solve  the   system  (3.4)  for  X  (P-,2^  ^"^^  ^    ^^12 
In  fact,  if  we  continue  this  procedure,  we  can  determine 
both  functions  along  the  entire  path.   Thus  we  have 
completely  determined  the  initial  data  for  these  functions, 
since  we  know  their  values  on  n  =  n,  by  reflection. 

Let  us  now  consider  the  interior  points  of  the  grid. 
In  general,  we  deteirmine  our  solution  point  by  point,  going 
from  left  to  right  in  each  row,  starting  with  the  bottom  row 
and  working  toward  the  top.   For  example. 

^11 '^12  ' ^22' ^13 '^^23' ^33 '^14 '^24 '^34 '^44 '^15'  *  *  * 

is  a  typical  ordering ;  note  that  we  stop  each  row  at  the 
diagonal,  since  the  remaining  points  are  known  by  symmetry 

(see  previous  section) . 

Suppose  now,  that  we  know  our  solution  at  P .  ,  .  and 
wish  to  find  it  at  P. ..   Presumably,  we  have   already  computed 
row    j-1  and  consequently  have  values  for  P.  .  ,  ,  as 
indicated  in  Figure  11. 

In  principle,  we  can  use  equations  (2 .14) -  (2 . 16)  in  the 
same  way  as  we  did  for  the  initial  data.   Unfortunately,  this 
requires  the  numerical  inversion  of  h(q),  which  is  very  time 
consuming.   Instead  we  find   u  and  v  by  solving  equations 

(2.17a,c)  numerically.   We  are  able  to  do  this,  because  we 
have  already  determined  their  values  on  the  initial  planes  by 
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the  previous  computation. 

To  be  more  specific,  we  wish  to  solve  the  following 
system: 

(3.5a)  V-  -  X^u  =  0 

(3.5b)  xj  +  X_Y^  =  0  , 

on  n  =  constant,  and 


(3.5c)  V  -  X  u  =0 

n    -  n 

(3.5d)  xjj  +  X_^yJ  =  0  , 

on  £  =  constant  (see  Section   2.4) ,  for  u(P. .) ,  v(P. .) , 
X^(P.  .)  ,  X-'-(P.  .)  and  Y-'-(P.  .)  . 

We  will  determine  u(P. .)  and  v(P. .),  and  therefore 
X^(P. .),  by  use  of  a  predictor  corrector  scheme.   Using 
equations  (3.5a,c),  we  construct  the  difference  equations 

(3.6a)   v(P..)  -  v(P._^^j)  =  X^(Pi_i,j)[u(P,.)  -  u(P,_,^j)]  , 

and 

(3.6b)   v(P..)  -  v(P.  .  ,)  =  X_(P    ,)[G(P..)  -  u(P   .   )]  , 

XJ  -L  /  J   J-  -^  f  J      -L        J-J  •*- f  J      -^ 

which  we  solve  for  the  unknowns   u(P. .)  and  v(P. .) .  Observe 
that  since  we  do  not  yet  know  X^(P.  .),  we  have  had  to  use 
X  (P.  ,  .)  and  X  (P.  .  ,)  which  implies  that  u  and  v 
approximate   u  and  v   to  first  orc^er  accuracy.  We  next  compute 
X^(P..),  based  on  u(P..)  and  v(P..),  and  use  them  in  (3.6) 
to  get  the  corrected  system 
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(3.7a)   v(P.  .)  -  v(P._^^.)  =  X^[u(P.  .)  -  u(P._^^.)]  , 


(3.7b)   v(P..)  -  v(P.  .  .)  =  X  [u(P..)  -  u(P.  .  .)]  , 


where 


and 


.   _  >.(p,^)  *  >.(p,,i.i) 


After  solving  the  system  (3.7)  for  u (P .  .)  and  v(P.  .) 
and  then  making  a  final  determination   of  X^(P. .),  we  can 
solve  the  difference  equations  which  result  from  equations 
(3.5b,d) , 


(3.8a)   X-'-(P.  .)-X-^(P.  ,  .)  + 
ij       1-1, J 


^-(Pij)^^-(Pi-l,j) 


[y-'-(p.  . )-¥■'•  (P.  .  .)] 
=  0 


and 


C^l 


(3.8b)   5i-^(P.  .)-X-^(P.  .  ,)  + 
13      i/D-1 


'X^(P.  .)+X^(P.  .  .) 
+   ij    +      1,3-1 


[y^(p..)-y1(p.^._,)] 

=  0, 


to  get  X""-  (P  .  . )  and  $■"■  (P  .  . )  . 
^       1]  1] 

We  state  for  once  and  for  all  that  this  is  the  method 

/si   ^1   ^2   '^2 
used  to  find  X  ,  Y  ,  X  ,  Y   and,  with  slight  modification, 

3      3 
X  and  Y  at  all  points  of  the  grid  for  which  they  are  defined. 

It  will  always  be  successful,  except  in  regions  where  the 


system  is  ill-conditioned,  i.e.  when 
1   X 


1   X 


=  X_^  -  X_  =  0  , 


q  =  c 


Thus  our  scheme  will  break  down  near  the  sonic   line, 
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4.    Calculating  the  Flow  in  the  Supersonic  Region 

In  this  section  we  explain  the  basic  reasoning  behind 
our  choice  of  supersonic  paths. 

The  essential  difficulty  with  the  calculation   in  the 
supersonic  region  is  that  although  we  want  our  integration 
paths  to  start  at  E,    ,g    ,^^,r)    ,t)    ,    and  n^  ,  in  keeping  with 
our  choice  of  initial  characteristics,  we  must  avoid  all 
points  of  {^,ri)  space  at  which 

since  our  difference  equations  will  be  ill-conditioned  there, 

and  the  procedure  will  break  down. 

To  begin  with,  we  consider  the  geometry  of  the  complex 

four-space  (?,ri)  (see  [30]  for  a  more  detailed  account). 

We  first  observe  that  there  is  a  two  dimensional  set  of 

points  which  we  shall  call  the  complex  sonic  surface,  whereon 

2    2^ 
q   -  c   =0 

(q  and  c  complex) .   This  surface  intersects  the  C  =  n  plane 
along  the  sonic  line,  whose  interior,  we  recall,  maps  into 
the  region  of  subsonic  flow  in  the  physical 

plar.e.       Observe  that  the  real  supersonic  regime  in  {?,ri) 
space,  though  attached  to  the  C  =  n  plane  along  a  portion  of 
the  sonic  line,  which  we  shall  call  the  real  sonic  segment, 
is  not  contained  in  that  plane,  since  for  q  >  c  ,  h(q)  is 
complex,  and  hence 
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s  =  h(q)e"^®  i-   h(q)e"^®  =  h(q)e^®  =  s*  , 
which  implies  that 

In  fact,  standard  theory  shows  ([  9  ],  pp.  526-535)  that  the 
characteristics  in  the  supersonic  regime  are  real,  and  all 
terminate  on  the  sonic  line. 

The  specific  task  at  hand,  then,  is  to  construct 
integration  paths  which  will  enable  us  to  compute  the 
solution  at  points  in  the  supersonic  region.   Put  another  way, 
we  want  to  choose  integration  paths  which  start  in  the 
subsonic  region  of  the  C  =  n  plane,  and  twist  around  the 
complex  sonic  surface  to  get  to  points  corresponding  to  the 
real  supersonic  region. 

To  do  this,   we  use  the  sonic  loci.   In  the  1=11, 
plane,  there  is  a  segment  of  the  (upper)  sonic  locus  v;hose 
points  are  connected  to  the  real  sonic  segment  of  the  C  =  n 
plane  by  5  characteristics.   Let  us  connect  this  segment  to 
the  initial  characteristics  by  a  path  P   (see  Figure  12A)  , 
choosing  a  net  of  points 


where  C-  /  i  >  k,  are  points  along  the  sonic  locus.  The  same 
situation  occurs  in  the  C  =  Cj^  plane,  where  we  choose  a  path 
Pt  (see  Figure  12B)  and  net 


'^l''^2'***''^i''^'i  +  l''^"i+2'***'^M  ' 
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subject  to  the  following  restrictions: 

(i)    we  require  that  for  every  ^-characteristic  passing 
through  P    and  n-characteristic  passing  through  P   ,  we 
have 

(3.9)  x_^(?,n)  i^   x_(c,n)  . 

Thus,  we  will  have  avoided  the  complex  sonic  surface, 
and  our  numerical  procedure  will  work  up  to  the  sonic  segment. 

(ii)   we  require  that  n  •  ,n  •  . -i  f  •  •  .  »1w  be  conjugate  points 
of  -k'^k+i ' • • • '?N  ^^   reverse  order,  i.e. 

(3.10)  nj  =  ?j^  ,  rij^i  =  fjj-l'  •••'  ^M  =  ^1  ' 

where  we  have  assumed  that  each  net  has  the  same  number  of 
points  on  the  sonic  locus,  or  rather,  that 

N-k+1  =  M-j+1  . 

We  can  best  understand  these  restrictions  by  looking  at  the 
resultant  grid  (see  Figure  13). 

We  first  remark  that  the  circled  diagonal  points  in 
the       diagram  are  in  the  n  =  C  plane,  by  relation  (3.10) 
and  further,  because  of  the  position  of  their  preimages  on 
the  sonic  loci,  that  they  are  points  on  the  real  sonic  segment. 
Second,   we  state  that  our  numerical  method  of  solution  (see 
previous  section)  works  here,  since  (3.9)  holds  everywhere, 
up  to  the  diagonal.  Of  course,   we  cannot  go  beyond  the 
diagonal,  since 
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there.   Observe,  also,  that  since  our  paths  are  asynunetric , 
we  no  longer  obtain  real  subsonic  points  on  the  main  diagonal. 

We  now  corr.e  to  the  entire  point  of  this  procedure, 
which  is  that  the  triangular  submatrix  of  points  beneath 
the  diagonal  (see  the  figure)  belongs  to   the  real  supersonic 
regime.   We  recall  that  the  surface  in  four-dimensional 
space  which  corresponds  to  the  real  supersonic  domain, 
is  attached  to  the  n  =  I  plane  along  the  sonic  line,  and 
further,  that  it  is  swept  out  by  characteristics  which 
terminate  at  the  real  sonic  segment.   Thus,  we  might  expect 
that  E,,,...,^      and  n  • /I  •,■./.•• /H..  together  sweep  out  the 
desired  region,  or,  put  another  way,   that  the  points  under 
discussion  yield  solution  values  in  the  supersonic  regime. 

We  note  that  the  preimage  of  the  real  sonic  segment 
is  best  located  by  observing  the  pattern  of  body  points  in 
the  subsonic  region  after  the  blade  has  begun  to  take  shape. 
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5.     Using  the  Existing  Method 

In  this  section,  we  discuss  attempts  at  using  the 
direct  analogue  of  the  airfoil  program  in  the  design  of 
compressor  and  turbine  blades. 

(A)    A  Compressor  Blade 

Dr.  D.  Horn's  first  attempt  was  relatively  successful. 
The  blade  in  Figure  8a,B  is  closed  and  of  a  reasonable  shape, 
though  it  has  a  high  gap-to-chord  ratio  and  may  suffer  from 
boundary  layer  separation.   Subsequent  blades  have  been 
designed,  and  are  undergoing  wind  tunnel  testing  now. 

It  should  also  be  noted  that   this  result  was  obtained 
by  using 


and 


s  =  1  -  2B5  -  C^ 


*        -      2 

s  =  1  -  2Bn  -  n 


for  equations  (2.16),  where  the  imaginary  part  of  B 
controls  the  lift  ([2  ]/  P-  69).  This  results,  of  course, 
in  a  5-plane  configuration  quite  similar  to  that  of  an 
airfoil. 

(B)    An  Unsuccessful  Turbine  Blade 

The  method  failed  in  the  turbine  case  (Fig.  14).  This 
was  not  totally  unexpected,  since  the  blade  was  of  such  high 
camber  and  large  thickness  as  to  be  unlike  anything  done  before. 
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The  choice  of  constants  in  equations  (2.16)  was,  here  and 
in  the  successful  configuration,  essentially  determined  by 
prescribing  the  graph  elements  (see  Section   4.4),  e.g. 
the  locations  of  the  nose,  the  branch  points,  etc. 

Although  the  placement  of  a  branch  point  at  (1,0) 
seemed  successful  in  magnifying  the  region  near  ?^  to  allow 
more  room,  for  the  placement  of  automation  paths  between  it 
and  the  sonic  locus,  we  were  still  unable  to  produce  a  closed 
blade.   It  appeared  that  the  lower  surface  had  an  unnecessarily 
high  curvature,  v/hich  caused  it  to  cross  the  upper  surface. 
Obviously,  this  case  was  physically  vmreasonable . 

Based  on  the  airfoil  experience  and  hundreds  of  runs, 
it  had  become  clear  that  the  method  failed  in  the  case  of 
extremely  thick  and  highly  camiered  blades. 
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4.  A  FURTHER  REFINEMENT  OF  THE  CHARACTERISTIC  TRANSFORMATION 

In  this  chapter,  we  return  to  equations  (2.16)  of 
Section  2.2,  which,  we  recall,  were  arrived  at  by  mapping 
flow  about  the  unit  circle  onto  its  interior.   Here,  we 
derive  a  new  expression  for  the  complex  velocity  by  instead 
mapping  the  flow  about  a  Joukowski  blade  onto  the  interior 
of  the  unit  circle.   Using  this  result  as  a  basis,  we  modify 
the  configuration,  or  equivalently ,  choose  new  constants  in 
(2.16)  and  obtain  a  turbine  blade. 

1.     The  Joukowski  Mapping 

This  section  shows  how  the  desired  physical  properties 
of  the  flow  (blade  thickness,  camber,  and  velocity  at 
infinity)  are  used  to  determine  a  specific  Joukowski  mapping. 

By  means  of  the  transformation  ([24],  pp.  194,195) 


(4.1)  2  =  i 


-.]■ 


a  pair  of  circles  tangent  at  ^  =  1  (see  Figure  15)   are 

mapped  into  the  z  plane  (see  Figure  17)   in  such  a  way  that 

the   circle     y^   =    ^^  =    ^   =   ''^0   "*"   ^^      ,    0   ;^  9    <    2Tr}    is  mapped 

onto  the  circular  arc   C   ,  and  circle   Yg=  (C:  C=CQ+be   ,  O<_0<  2-n} 

is  mapped  onto  the  Joukowski  profile  C  . 

For  a  flow  about  C„  with  circulation  T   and  complex 

B 

velocity  U  at  infinity,  the  com.plex  potential  is  given  by 

,2      .„      C-Co 
(4.2)         X(c)  =  U(c  -  Cq)  +  U(^)  +  j^  In  [^^—]     . 
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Using  (4.1)  and  (4.2),  we  get  the  complex  velocity. 


w(r)  =  ^  =  ^  ^  = 
^^^'         dz    di;  dz 


Ub^        .,        .,2 


U 2  -^ 


(Z-^q)  2tt(C-;q)J  (C  -1) 


2C 


[2ttU(C-Cq)^+  ir(;-CQ)  -  27TUb^]C^ 
^  2 2 * 

In  addition,  the  Kutta  condition  for  a  cusped  tail 
requires  that  c  =  1  be  a  root  of  the  numerator,  so  that 

2;^[U(C-Cn)  +  U(l  -  l^)] 

(4.3)  w(C)  = 5 —   , 

(C+1)  U-^q) 

and 

(4.4)  r    =      4TT[uQeo    +    (1    -   ctQ)vQ]     , 

where  we  have  set  Cq  =  Cq  "•■  ^^ q   ^^^   "^  ~  "o  ~  "^^0' 

The  flow  will  be  completely  determined  by  equations  (4.1)- 
(4.4)  once  we  have  properly  chosen  the  as  yet  undetermined 
constants.   In  practice  we  specify  the  velocity  at  infinity  U, 
the  velocity  angle  at  the  tail  G   ,  and  a  thickness  parameter  t, 
which  are  then  used  to  determine  t. ,  a,  Cq  t    and  b.   Because 
of  the   symmetry  of  the  circular  arc,  9   is  essentially  the 
blade  camber.   In  the  following,  we  assiome  that  U  and  6   are 
given. 

The  complex  velocity  at  the  tail,  w_,  =  u„  -  iv„  is  given 

by 

2  Re  (D(1-Cq)) 
w_,  =  w(1)  = 


so  that 

tan 


(I-Cq)" 


v^   -  Im  [(I-Cq)^]       2BQ(aQ-  1) 


■^   ""t     Re  [(1-Cq)^]       (1-0^)^-60 
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Solving  for  B-,  we  get 


Bo  =  (l-^o) 


cos  0^  ±  1 


sin 


These  two  equations  define  a  pair  of  straight  lines 
in  the  ^-plane  along  which  c,^   must  be  chosen  to  get  the 
desired  6  .   In  our  configuration  we  have  assumed  6   <  0, 
a.  <  0,  and  6f^  >  0  which  requires  the   choice  of  the  minus 
sign  in  the  above  (see  line  I   in  Figure  16 ) •   The  location 
of  C«  on  this  line  will  be  determined  by  the  thickness 
parameter  t. 

Ktxt,  we  determine   t„  and  a,  the  center  and  radius, 
respectively,  of  ■)■.    Since  y,   ^^  tangent  to  Yq  at  C  =  1, 
and  since  y.    also  passes  through  c  =  -1,  it  is  clear  that 
t-  must  be  at  the  intersection  of  line  £  and  the  imaginary 
axis.   Hence, 


and 


^0  = 


a  =  Il-t, 


cos  e,p  -  1 


sin 


/2(l-cos  e^) 


sm 


We  now  define  the  thickness  parameter  t  as 
t  =  b  -  a  >  0  , 
which  is  the  distance  of  c,^    from  t-  along  I.    If  t=0,  i.e.  we 
have  the  case  of  zero  thickness,  y      =  Yj,  and  the  profile 
coincides  with  the  circular  arc.  As  t  increases,  the  body 
will  thicken  about  the  arc. 

Using  the  above  equations,  we  can  now  generate 
Joukov/ski  profiles  by  specifying  U,  G„  and  t. 
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2 .    Mapping  the  Flow  onto  the  Interior  of  the  Unit  Circle 

In  this  section  we  map  the  flow  in  the  C-plane  onto 
the  interior  of  the  unit  circle  in  the  5-plane,  thus  determin- 
ing the  functions   z(^)  and  w(^). 

This  is  accomplished  by  means  of  the  transformation 

(4.5)  5  =  ^^, 
or,  solving  for  t,, 

(4.6)  ^  =  ^^  +  I  . 

Substituting  (4.6)  in  equations  (4.1)  and  (4.3),  we  get 


.„-^*    ' 


and 

2(b+CCn)^KU(l-Cn)+Ub] 
(4.8)        w(0  =  


[^(I+^q)  +  b]b^ 


Clearly  then,  w  is  a  rational  function  of  ?  which,  in 
addition  to  its  branch  points,  has 


(1)    a  pole  at 


(2)    a  zero  at 


5=  -" 


1*^0  ' 


5  =  r^ 


which  corresponds  to  the  nose  of  the  blade,  and 

(3)    a  double  zero  at 

5  =  -^ 


^0 
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The  new  elements  introduced  here  are  the  pole  and  the 
double  zero.   The  Joukowski   transformation  suggests  that 
the  location  of  the  pole  is  primarily  determined  by  the 
thickness  of  the  blade,  whereas  the  position  of  the  double 
zero  is  primarily  determined  by  the  camber  (see  Section  3) . 
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3.     Modifying  the  Configuration 

In  this  section  we  see  how  the  physical  characteristics 
of  the  Joukowski  blade  affect  the  configuration  in  the 
5-plane,  viz.   the  location  of  the  pole,  double  zero,  and 
branch  points  of  w. 

The  three  parameters  of  interest  are:  (a)  the  thickness 
of  the  blade,  (b)  its  camber,  and  (c)  the  angle  of  the  flow 
at  infinity,  which  respectively  correspond  to  t,  6   and  U 
of  Section  1.   Note  that,  in  the  following  figures,  an 
additional  rotation  has  been  introduced  so  that  the  nose  is 
always  located  at  -1. 

(a)  Blade  thickness   (Figure   18) 

Here  we  vary  t  while  keeping  G   and  U  constant. 

As  can  be  seen  in  the  figure,  the  blade  becomes 
thicker  as   t   increases,  and  the  pole  moves  away  from  the 
nose.   The  double  zero  is  off  the  page  and  relatively 
unimportant. 

(b)  Blade  camber     (Figure   19) 

Here  we  vary  Q      while  keeping  t  and  U  constant. 

As  6   increases  negatively,  the  blade  camber 
increases,  the  pole  approaches  the  nose,  and  the  double  zero 
comes  into  view. 
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(c)    Angle  of  the  flow  at  infinity   (Figure  20) 

We  fix  t,  e   ,  and  |u|  =1,  while  varying  the 

argument  of  U. 

This  in  no  way  changes  the  shape  or  size  of  the  body 

in  the  z-plane,  but  drematically  affects  the  imaginary 

coordinate  of  the  pole,  as  well  as  the  relative  positions 

of  the  branch  points. 

Our  procedure  then  is  to  use  the  Joukowski  blade  to  get 

a  preliminary   estimate  for  the  configuration,  and  then  to 

vary  the  locations  in  accordance  with  our  design  goals. 
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4 .    Conclusion;  A  Turbine  Blade 

In  Figure  21,  we  exhibit  a  turbine  blade.   Observe  that, 
as  a  result  of  our  particular  choices  for  the  constants  in 
equations  (2.16),  the  configuration,  especially  the  sonic 
line  (or  locus) ,  is  quite  different.   Furthermore,  it  happened 
that,  rather  than  choose  the  C.  directly,  it  v:5s  easier  to 
determine  them  by  prescribing  certain  basic  elements  of  the 
transformation . 

In  this  case,  we  assumed  the  equations  to  be  of  the 
fourth  degree,  which  meant  that  we  could  satisfy  five  condi- 
tions in  the  C-plane  (Figure  21B) . 

(1)    We  put  the  leading  edge  stagnation  point,  the  nose,  at 
-1+Oi,  i.e.   we  required 

(4.9a)  C.4  =  0 


where 


Cn  =  -i^oi' 


(2)    We  placed  a  branch  point  at  ^^p  =  1+Oi  to  make  room 
between  the  body  and  the  sonic  line,  i.e.  we  required 


(4.9b)  ^  (C,C^) !      =  0 

•HP 


dC   X   1^=, 


(3)    Further,  we  specified  the  velocity  at  that  point  (a  Mach 
number  of  .9  and  a  flow  angle  of  148.93  degrees),  which  led 
to  a  second  equation. 
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(4.9c)  C.^Jp  =  h(q3p)  e    ^^  . 

(4)    Finally,  we  required  a  double  zero,  as  dictated  by  the 
Joukowski  transformation,  which  we  placed  at 

K^  =   -.35  +  1.4i  . 

This  led  to  the  two  equations, 

(4.9d)  C^CJ  =  0  , 

and 

=   0  . 


(4.9e)  ^  (C^?^) 


^  =  ^D 


The  last  requirement  is,  of  course,  the  essence  of  this  work. 
Note  that  ^  is  in  the  same  general  area  as  in  the  incompres- 
sible case  (Figure  19D) . 

We  remark  that,  having  determined  the  C.  by  the  five 

equations  (4.9),   the  locations  of  E,      and  £„  were  found  as  a 

A      B 

consequence  of  the  prescribed  inlet  and  exit  velocities. 

Thus,  though  the  pressure  distribution  along  the  blade 
indicates  probable  boundary  layer  separation,  it  is  clear  that 
the   method  works,  i.e.  we  can  design  highly  cambered 
blades  of  large  thickness.   It  is   hoped  that  a  physically 
practical  blade  will  be  designed  and  tested  in  the   near  future, 
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This  report  was  prepared  as  an  account  of 
Government  sponsored  work.   Neither  the 
United  States,  nor  the  Administration, 
nor  any  person  acting  on  behalf  of  the 
Administration : 

A.  Makes  any  warranty  or  representation, 
express  or  implied,  with  respect  to  the 
accuracy,  completeness,  or  usefulness  of 
the  information  contained  in  this  report, 
or  that  the  use  of  any  information, 
apparatus,  method,  or  process  disclosed 
in  this  report  may  not  infringe  privately 
owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to 
the  use  of,  or  for  damages  resulting  from 
the  use  of  any  information,  apparatus, 
method,  or  process  disclosed  in  this 
report . 

As  used  in  the  above,  "person  acting  on  behalf 
of  the  Administration"  Includes  any  employee 
or  contractor  of  the  Administration,  or 
employee  of  such  contractor,  to  the  extent 
that  such  employee  or  contractor  of  the 
Administration,  or  employee  of  such  contractor 
prepares,  disseminates,  or  provides  access  to, 
any  Information  pursuant  to  his  employment  or 
contract  with  the  Administration,  or  his 
employment  with  such  contractor. 
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